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Abstract 

We study the effective dynamics of a mechanical particle coupled to 
a wave field and subject to the slowly varying potential V(eq) with 
e small. To lowest order in e the motion of the particle is governed 
by an effective Hamiltonian. In the next order one obtains "dissipa- 
tive" terms which describe the radiation reaction. We establish that 
this dissipative dynamics has a center manifold which is repulsive in 
the normal direction and which is global, in the sense that for given 
data and sufficiently small e the solution stays on the center manifold 
forever. We prove that the solution of the full system is well approxi- 
mated by the effective dissipative dynamics on its center manifold. 



1 Introduction 



At the beginning of this century, in the context of the Maxwell-Lorentz equa- 
tions, radiation reaction was one of the most outstanding problems in the- 
oretical physics. It was left sort of unfinished when theoreticians turned to 
quantum electrodynamics. In this paper we study radiation reaction in the 
mathematically somewhat more accessible case of a scalar wave field. We 
believe that our results provide good indications on the effective dynamics 
for a charge coupled to the Maxwell field ^2 . 



To explain in more detail the physical context we have to set up the model 
first. We consider a particle, position q(t) G IR 3 and momentum p(t) G IR 3 , 
with "charge" distribution p of total charge 

e = J d 3 xp(x) ^ 0. 

We require that p is smooth, radial, and supported in a ball of radius R p , 

p G C °°(IR 3 ) , p(x) = p r (\x\) , p(x) = for |ac| > R p . (C) 

The particle is coupled to the scalar wave field 4>(x, t) with the canonically 
conjugate momentum field n(x, t), x G IR 3 . In addition the particle is subject 
to an external potential, V, whose properties will be listed below. We assume 
that the potential is slowly varying on the scale of the charge distribution, i.e., 
on the scale set by R p . Formally we introduce the dimensionless parameter 
£, £ « 1, and consider the scale of potentials V(eq), e — ► 0. The equations 
of motion for the coupled system are 

<p(x, t) = n(x, t), 7r(x, t) = Acj)(x, t) — p{x — q(t)), 

v(t) r (i- 1 ) 

q(t) = , F U , p(t) = -eW(eg(t)) + / d 3 x(p(x,t)\7p(x - q(t)) . 
'l+p(t) 2 J 



The dynamics governed by (|1 . 1|) has three distinct time scales, well- 
separated as e — > 0. On the microscopic time scale, t = 0(1), the particle 
moves along an essentially straight line and the field adjusts itself stationar- 
ily On a time scale that we call the macroscopic scale, the particle 

feels the potential and responds to it with an effective kinetic energy which 
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incorporates the coupling to the field. This scale was studied in H . The par- 
ticle looses energy through radiation at a rate roughly proportional to q(t) 2 . 
Thus on the macroscopic time scale, friction through radiation is of order e. 
To resolve such an effect we have to go to even longer times or to look with 
higher precision. The friction time scale is the subject of our paper. 

The dynamics of ( |1 . 1| ) is of Hamiltonian form. We need a few facts in the 
case the external potential vanishes, V = 0. Then (|1 . 1| ) has the energy 



H (<P,n,q,p) = (l+ p 2 ) 1/2 + ^ / rf 3 x(|vr(x)| 2 +|V0(x)| 2 )+ J d 3 x <P(x)p(x - q) 
and the conserved total momentum 

V{4>, 7T, q, p) = p + J d 3 x(J)(x)\7-K(x). 

The minimum of Tio, at fixed V, is attained at 

S q ,v = (<f>v(x -q),rr v (x- q),q,p v ) 



;i.2) 



where v G V = {v : \v\ < 1}, p v = v/yl 



-v ■ V0„, and <p v (k) 



—p{k)/[k 2 — {y ■ k) 2 }; the hat denotes Fourier transform. We call S qtV the 
soliton centered at q, v. It has the normalized energy 



S s (v) = Tto(S qjV ) — Ho(S qj0 ) 

= (1 -v 2 Y 1/2 - 1 + 3m e 

and the total momentum 

V s (v) = P(S q , v ) 



2 1 
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;i.3) 



Here m e = | / <i 3 fc |p(/c)| 2 /c~ 2 is the mass of the particle due to the coupling 
to the field. We note that because of the Hamiltonian structure we have the 
identity v(dV s /dv) = (d£ s /dv). It is shown in that the map v V s (v) is 
invertible from V to IR 3 , with inverse P \— > v(P), and thus E(P) := £ s (v(P)) 
is well defined. 
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Taking S q%v as initial conditions for ( |1.1| ) with V = we obtain a solution 
travelling at constant velocity v, 

S q , v (t) = (<t>v(x -q- vt),7T v (x - q- vt),q + vt,p v ) , v G V. 

Let us call {S q>v : q G JR 3 ,v G V} the six-dimensional soliton manifold, S. 
Thus, for V = 0, if we start initially on S the solution remains on S and 
moves along the straight line t i— > g° + f /f:. In fact, if we start close to S, then 
S is approached asymptotically, 0. When the particle is subject to a slowly 
varying external potential, then the rough picture is that the solution will 
remain close to S in the course of time. For simplicity we assume throughout 
that the initial datum for ( |1 . 1|) lies exactly on S, i.e., 

(0(O),vr(O),g(O),p(O)) = ^o il ,o, (1.4) 

possible generalizations being discussed below. 



At this point it is instructive to transform fll.lj) to the macroscopic space- 
time scale in such a way that the field energy remains constant. Then the 
macroscopic variables, denoted by a ', are 

t = £ -H', q = £ - 1 q', x = e~ 1 x', q(t)=e- 1 q'(t'), 
and (j)(x,t) = y/i(p'(x', t'). 

We also set 

p £ (x) = e^pie^x) . 

In particular, p £ (x) = for \x\ > eR p and / d 3 xp £ (x) = J d 3 xp(x). With this 
convention, omitting the primes and indicating explicitly the ^-dependence 
of q'{t'), we arrive at 

4>(x,t) = A(j)(x,t)-^p £ (x-q e (t)), (1.5) 
q £ (t) = v £ (t), 

m Q {y £ {t))v £ {t) = -VV(q(t)) + y/e f d 3 x <j>(x, t)Vp E {x - q £ {t)). 



Here mo(v) is the 3x3 matrix defined through mo(v)v = 7?) + 7 s (t> • v)v 
with 7(f) = 1 /yl — v 2 . Rather than momenta as in (|1 . 1| ) , we use velocities 
which turns out to be more convenient in our context. The initial soliton 
( |1.4|) transforms to 

S e q o, v o = (fAx-q°),< (x-q ),q°,v ), (1.6) 
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where (j)%{k) = —^/ep{ek)/[k 2 — (v ■ k) 2 } and 7r £ = — v ■ S7(j) £ v . Thus, on the 
macroscopic scale, the total charge is ^fe J d 3 xp(x), whereas 



±e J d 3 k\p £ (k)\ 2 k- 2 
is independent of e. Eqs. ( |1.5[) are again of Hamiltonian form. The energy 
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U mac {^^v) = l(v) + V(q) + -J d 3 x(\rc(x)\ 2 + |V0(x)| 2 ) 

J d 3 xcj){x)p £ {x - q) (1.7) 

is conserved under ( |1.5| ). It is bounded from below, as 7i mac (0, 7r, g, t>) > 
V(g) — 3m e independently of e. 

There is another, very instructive way to think about the initial value 
problem (|1.5|) , ( |1.6|) . We prescribe initial data at t — —r, r > 0, which have 
finite energy and some smoothness. We refer to for the precise conditions. 
We solve (|1.5|) for V = up to time i = 0. Then in the limit r — > oo the data 
at t = are exactly of the form (|1.6| ). For t > the external forces are acting. 
Clearly this causes some mismatch, which is reflected by a non-smoothness 
of the fields (0, 7r) at the light cone {x : \x\ = t, t > 0} in the limit e — > 0. 

Under suitable assumptions on V and for |p| L2 sufficiently small we proved 
in [§ that 

|<f(t)|<^<l, \q £ (t)\<C, and | g' £ (t)| < C (1.8) 
uniformly in e and t G E, and that the limit 

lim q £ {t) = r(t) (1.9) 

exists. Here r(t) is the solution of Hamilton's equations of motion with the 
effective Hamiltonian E(p) + V(q), cf. the definition of E(jp) below ( |1.3| ), 
which in terms of velocities read 

r = u, m{u)u = — VV(r), (1-10) 

with initial data r(0) = q°, u(0) = v°. Here m(u) = mo(ti) + rrif(u), where 
rrif(u) is the additional "mass" due to the coupling to the field defined by 

mf(u)u = 3m e (tp(\u\)u + |w| _1 (/? / (|w|)(w • ii)uj (1-H) 
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as a 3x3 matrix, where y(H) is the function appearing in the square brackets 
of Eq. (|1.3| ) . Note that the energy 

H(r,u) =£ s (u) + V(r) (1.12) 

is conserved by the solutions to (|1.10Q . 

With this background information let us return to the radiation reaction 



as discussed by Abraham, Lorentz, Schott, and Dirac, cf. [|14| for an excellent 
account. Of course, these theoretical physicists were interested in the elec- 
trodynamics of moving charges. We take here the liberty to transcribe their 
arguments to the case of a scalar wave equation. For the sake of discussion 
we reintroduce the bare mass m and state the equations for small velocities 
only. In our proof below, however, we will handle all v 6 V. 

At the beginning of this century the hope was to define a structureless 
elementary charge through a point charge limit. For this program, one had 
to model the charge distribution phenomenologically with the understanding 



that finer details should become irrelevant in the limit. In ( |1.5|) we adopted 



the Abraham model of a rigid charge distribution. The point charge limit 
then corresponds to taking in (|1.5f) the charge distribution p e instead of \fsp e - 
With this choice e = J d 3 xp £ (x) = J d 3 xp(x), whereas the electromagnetic 
mass equals | / d 3 k\p E (k)\ 2 k~ 2 = e~ 1 m e . A formal Taylor expansion leads to 
the effective equation of motion 

m r = — VV(r) — e~ x m e r + ae 2 r, (1-13) 



valid for small velocities r, with some constant a > 0. Eq. ( |1.13| ) is the non 



relativistic limit of the Lorentz-Dirac equation, [jnj. The standard argument, 
reproduced in many textbooks, e.g. ||, (with the notable exception of Lan- 
dau and Lifshitz ||) is to lump mo and e~ 1 m e together and to take the limits 
e —>■ and mo — > — oo at constant m + e~ x m e = m exp , the experimentally 
observed mass of the particle. Then ( |1.13| ) reads as 

m exp f = -W(r) + ae 2 r . (1-14) 

Since this equation is of third order, one needs besides q°, v° also u(0) as 
initial condition which has to be extracted somehow from the initial data 
of the full system. Even worse, (|1.14 ) has solutions which are exponentially 



unbounded in time, the famous run-away solutions. Thus one needs an addi- 
tional criterion to single out the solutions of physical relevance. Dirac ffl, and 
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later Haag [f2], argued that physical solutions have to satisfy the asymptotic 
condition 

limr(t)=0, (1.15) 

t — *oo 

as a substitute for the missing initial condition f(0). The validity of the 
asymptotic condition has been checked only in trivial cases; see [plj. For 
general V one should expect the solutions to Q1.13 ) to be chaotic. Physical 



and unphysical solutions might be badly mixed up. On a more practical 
level, the physical solutions are unstable and therefore difficult to compute 
numerically. To put it in the words of W. Thirring jn|: u ... (\l. 14) has not 



only crazy solutions and there are attempts to separate sense from nonsense 
through special initial conditions. But one hopes that the true solution to the 
problem will look differently and that the nature of the equations of motion 
is not so highly unstable that the act of balance can be achieved only through 
a stroke of good fortune in the initial conditions." 

This is indeed the case, as we are going to show in this paper, and our 



resolution requires just a little twist. If according to ( |1.5| ) we adopt the 
macroscopic time scale, then ( |1.13| ) reads 

(m + m e )r = — W(r) + eae 2 r (1-16) 

which just reflects that radiation reaction is a small correction to the Hamil- 
tonian motion. The bare mass tuq should be kept strictly positive. Otherwise, 
TL mac from ( |1.7|) is not bounded from below and ( pT5|) has solutions increas- 
ing exponentially in time, a phenomenon completely unrelated to run-away 
solutions, however. 

In (|1.16| ) the highest derivative appears with a small prefactor. Such 



differential equations are studied in geometric singular perturbation the- 
ory. From there we know that ( |1.16| ) has a six- dimensional invariant cen- 
ter manifold X £ , which is only 0(e) away from the Hamiltonian manifold 
1q = {(q,q,q) : (m + m e )q = — W(q)}. For initial conditions slightly off 
I e the solution moves away from X £ exponentially fast. On X £ , q is bounded 
away from 1, q is bounded, and the motion is governed by an effective second 
order equation, cf. Eq. (fO|) below, which gives precisely the physical solu- 
tions. To establish such a result we have to prove that the solution to ( |1.5| ) 
stays indeed close to I £ . 

In our paper we carry out this program, essentially under the same con- 
ditions as in ||, namely a sufficiently differentiable V and \p\ L 2 small. Our 
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main additional estimate is 



\v (t)\<C (1.17) 

uniformly in e and t G IR. Thereby we can bound one further order in the 
rigorous Taylor expansion and obtain, setting q £ = v £ , 

m(v £ )v £ = -W(<f) + ea(v £ )v £ + eb(v £ ,v £ ) +e 2 f £ (t), t > et u (1.18) 

with |/ e (t)| < C and coefficient functions a, b that will be defined below. 
Clearly (|1.18|) should be compared with 

r = u , m{u)u = — VV (r) + ea{u)u + sb(u, u) . (1-19) 

Our crucial observation is that the condition \u(t)\ < const. < 1 for all t holds 
only on the center manifold X e . Thus the a priori estimate |<7 e (t)| < v < 1, 
see ( P-.8j) , together with the initial conditions r(0) = q°, u(0) = v°, uniquely 
singles out that solution of (|1.19|) which is to be compared with the true 
solution. 

Since on the error term f £ (t) in Q1.18D we only know that it is uniformly 



bounded, the difference \q £ (t) — r(t)\, with r(t) having initial conditions on 
X e , can be bounded at best as ee ct . Thus on the time scale t = 0{1) we seem 
to be back to the result ( |1.9| ) already proved in 0. To distinguish, from 
this point of view, between ( |1.19| ) and ( |1.10| ) we would have to control the 
difference with a precision of order e 2 . At present we do not know whether 
this is possible, but nevertheless we can prove the weaker statement 

\H(q £ (t),v £ (t)) -H(r(t),u(t))\ < const, e 2 , (1.20) 

where H is the energy from (1.12). Thus on a surface of constant energy 



the difference \q £ (t) — r(t)\ could be of order e, whereas along ViJ it must 
be of order e 2 . In addition to ( 1.20| ) it may also be shown that in fact 



\q £ {t) — r(t)\ ~ e 3 on the short time scale t = 0(e), a result that is quite 
natural from the viewpoint of singularly perturbed ODEs. On the original 
time scale of ( |1 . 1| ) this amounts at least to an estimate with precision e 2 over 
time intervals of length 0(1), a result that could not have been obtained 
from the bounds in 01. 
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2 Main results 



We give some more details and state our main results precisely. First we 
have to establish the bound (|1.17|) . 

Lemma 2.1 For \p\ L 2 sufficiently small we have 

sup | v £ (t)\<C 
ten 

for every solution of ( \l.dj ) which starts on the soliton manifold S. Both the 
constant C and the bound for \p\ L 2 depend only on the initial data. 

The bound of Lemma 2.1 may be used to Taylor expand the self- force 

F'(t) = Ve J d 3 x<f)(x,t)Vp £ (x-q £ (t)) (2.1) 

in ( |1.5| ) as 

F s e (t) = -m f (v £ (t))v £ (t)+ea(v £ (t))v £ (t)+eb(v £ (t),v £ (t)) + 0(e 2 ) , t>et u 

(2.2) 

which together with the second equation in (|1.5| ) yields ( |1.18|) . Here m/ is 



defined in ( |1 . 1 1| ) , and t\ = 2R P /(1 — v) is the microscopic time the wave 



equation needs to forget its data because of the compact support of p and 
the velocity bound, cf. assumption (C) and Q1.8Q. The coefficient functions 
are given by 

a{v)v = (e 2 /247r)(i)- V„)V„7 2 = (e 2 /127r) [ 7 4 i) + 4 7 6 (w • i))w] , (2.3) 
b(v,v) = (e 2 /32n)(v- V,) 2 V„ 7 2 

= {e 2 /4:7r)[2-f 6 (v-v)v + 7 6 v 2 v + 6"f 8 {v-v) 2 v}, (2.4) 



v,v e IR , with 7 = 1/Vl - v 2 , \v\ < 1. 

Next we explain the existence and the role of the center-like manifolds 
I £ in greater detail. We refer to JTl, for further background on geometric 



singular perturbation theory. To rewrite ( |1.19| ) as a singular perturbation 
problem, let 

x = (r,u) e H 3 x V, y = uelR 3 , f(x, y) = (x 2 , y) G V x IR 3 , and 
g(x,y,e) = a(x 2 )~ 1 [m(x 2 )y + VV(zi) -eb(x 2 ,y)) . 
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Then ( 



1.19| ) reads as 



x = f(x,y), ey = g(x,y,£). (2.5) 



We intend to apply the results from [Til] to (|2.5|) in order to find a center-like 



manifold for the perturbed problem near the corresponding manifold for the 
(e = 0)-problem. With h(x) = — m(x 2 ) _1 W(xi), let 

2b = {(x,y) ■ g(x,y, 0) = 0} = {(r,u,it) : m(u)tt = -W(r)} 

= {(i,/i(x)):x6E 3 xV} (2.6) 



be this invariant manifold for ( |2.5| ) with e — 0. The flow on Xo is governed 
by the equation x = f(x,h(x)), or stated differently, m(f)r = — VV(r), the 
familiar Hamiltonian flow. 

To see that Xo is perturbed to some X e with e small, we have to modify the 
functions a(u), m(u), and b(u,u) for \u\ close to one due to the singularity 
at \u\ = 1. This will cause no problems later on, since we already have the 
a priori bound |f e (t)| < v < 1 for the velocity of the true system. In (|4.4j) 
below, we will fix a small 5 = 5(v) > satisfying some estimates; 5 depends 
only on bounds for the initial data, since v does so. Let 

/d_j = IR 3 x {u E M 3 : \u\ < 1 - 5}, 

We continue a(tt), m(u), and b(u,u) with their values at |it| = 1 — 8 to the 
missing infinite strip 1 — 5 < \u\ < 1. Then the basic assumptions (/), 
(II) from p. 45] are satisfied, since Xo is also what is called normally 
hyperbolic, i.e. repulsive in the direction normal to Xo at an e-independent 
rate, see Lemma JO] below. Hence we find e = e {6) > and a C 1 -function 
h(x,e) = h £ (x) : IR 3 x Vx]0,e ] —> H 3 such that for e < e Q , 

I e = {(x,h E {x)) : x G IR 3 x V} 



is forward invariant for the flow ( 1.19 ) with the modified functions a, m, b 



Since the modified equation agrees with (|1.19|) in the interior of /C-^, we 
conclude that X £ is locally invariant for the flow ( |1.19|) , i.e. the solution of 



the modified equation is the solution to the original equation as long as it 
does not reach the boundary set {(x,h e (x)) = (r, m, h £ (r,u)) : \u\ = 1 — 5}. 
The flow for e = is then perturbed to x = f(x, h E (x)) for e < e . 
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We will show in Theorem below that for e G]0,£i], with e\ > 
sufficiently small, all solutions of ( |1.19|) starting at points (r, u, h s (r,u)) G T £ 



with \u\ < v, will indeed stay away from the boundary {(r, u, h £ (r, u)) : \u\ = 
1 —5} for all future times. In addition, W(r(t)) — > and r(t) — > as t — > oo, 
which is just the asymptotic condition (|1.15|) postulated by Dirac and Haag. 
If the potential is sufficiently confining, then the solution trajectory on I e 
not only approaches the set of critical points for V in the long-time limit, but 
it converges to some definite critical point. Moreover, we will show that for 
all solutions on the center manifold, u(t) and u(t) are bounded, and u(t) is 
bounded away from 1, uniformly in e and t. Conversely, every such solution 
to (|1.19|) has to lie on X £ . Thus X £ indeed characterizes the physical solutions. 



To summarize, we have established now the existence of a center manifold 
X e with a well-defined (semi-) flow on it that gives a unique solution to ( |1.19|) 
for initial velocities bounded by v. 

For the potential V G C 3 (IR 3 ) we assume that it is bounded in the sense 
inf g6R 3 V(q) > — oo and 

sup (\V(q)\ + \VV(q)\ + \VVV(q)\ + |VVW(g)|) < oo . (U) 

The method works equally well for V G C 3 (IR 3 ) which is confining, i.e., 

V(q) — > oo as \q\ — > oo , (U') 
as will be made more precise in Section fl, cf. Theorem PO. 



Our main result is the following 

Theorem 2.2 Assume (U) or (U') for the potential, and let the initial data 
(4>°(x), 7r°(x), q°, v °) for (\l.dj) be given by (\l.dj) . Let \p\ L2 and e < e\ be 
sufficiently small, and introduce the center manifolds I s for the comparison 
dynamics fli.ify ) as explained above. At time et\ = e2R p /(l — v) we match 
the initial values, r(et\) = q £ (eti), u(eti) = v £ (eti), for the motion on the 
center manifold, i.e., the initial data for the comparison dynamics are 

(q'ieh), v £ (eh), /i £ (<f (eti), v'ieh))) G 1 £ . 

Then for every r > there exists c(r) > such that for all t G [eti,eti + r] 

\q £ (t)-r(t)\<c(r)e, \v £ {t) - u{t)\ < c(r)e , and \v £ {t) - u(t)\ < c(r)e. 

(2.7) 
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In addition we have the bound 



\H{q^t)^{t))~H{r{t)Mt))\<c{r)e 2 . (2.8) 

Remarks 2.3 (i) As already mentioned at the end of the introduction, we 
can also show 

\q £ (t) - r(t)\ < c(t)s 3 and \v £ (t) - u(t)\ < c(r)e 2 (2.9) 

for t G [et\,et\ + et\, i.e., t = 0(s), cf. Proposition [57J]. 

(ii) The construction of the center manifolds and the upper bound for \p\ L 2 
rely only on bounds for the data, but not on properties of a particularly 
chosen solution. Our main technical assumption is a sufficiently small \p\ L 2 
which is presumably not necessary. 

(iii) In |J we did not require the true solution to start on the soliton manifold, 
but instead to start close to it. We refer the criterion [||, Thm. 2.6] for an 
"adiabatic" family of solutions. The same generality could be achieved in the 
present context, using an appropriately modified version of 0, Thm. 2.6]. 
In Section |8| we derive the relevant estimates, in particular (|8.8| ), in full 
generality containing a non-zero initial difference Z(0). The corresponding 



generalization of Theorem [272] is then straightforward. However, since we did 
not want to obscure our main achievement through technicalities, we decided 
to elaborate here the more accessible case of a trajectory starting right on 
the soliton manifold. In the same spirit we do not consider arbitrary time 
intervals of length r, but only the particular [eti,eti + r]. 

(iv) The existence of solutions to (|1 . 1|) is discussed in Lemma 2.2]. For 
every initial value Y° = (<f)°(x), vr°(x), q°,p°) G £ we find a unique (weak) 
solution Y(-) G C(TR, £) such that Y(0) = Y°. Here the state space is 
£ = .D^^IR^eL^lR 3 )©^ 3 ©^ 3 [where D 1 - 2 (B 3 ) = {<p G L 6 (IR 3 ) : \V(f>\ G 
L 2 (1R 3 )}] with norm \Y\ £ = \V(j>\ L2 + \n\ L2 + \q\ + \p\. 

Having such fairly precise information on the particle trajectory we can 



also determine the adiabatic limit e — > of the fields (0, tt) in (1^5) through 
the solution of the inhomogeneous wave equation. We generate the ini- 
tial data as explained in the introduction. On the level of the compari- 
son dynamics this means to extend r(t) and u(t) to negative times t < 
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by r(t) = q° + tv° resp. u(t) = v°. Let the retarded time t Tet , depend- 
ing on x and t, be the unique solution of t Tet = t — \x — r(t ret )|, and let 
n(x,t) = (x - r(t Iet ))/\x - r(t Iet )\. 

Theorem 2.4 Under the conditions of Theorem and for the fields ((f), ir) 
from we have for x ^ r(t) the pointwise limits 

lim —=(f)(x,t) = . —. — - (l-n(x,t) ■ u(t Tet )Y l (2.10) 



e^o ' 47rja; — r(tret)| 

and, except for the light cone {x : \x\ = t > 0}, 

lim —= ir(x, t) 



■— — — t (l - n(x,t) ■ u(t rct )) h(x,t) ■ u(t ret ) 

Atc\x - r(t Ict )\ v > 

-— I (1 - n(x,t) ■ u(t TCt )) [n(x,t) ■ u(t Tet ) - u(t 

4:1t\x - r(t rct )\ 2 v ' 



vet) ) 



(2.11) 



The paper is organized as follows. Since the proof of Lemma ^TT] is rather 
technical, we moved it to an appendix, Section |8|. The derivation of the 
representation ( |2.2| ) of the self-force term is the contents of Section [| In 
Section f| we give supplementary remarks on the behaviour of solutions on 
the center manifold, whereas in Section |^ we carry out the proofs of Theorem 
|2.2| and Proposition |5.1| . Section ^] contains the proof of Theorem [2.4| , and 
finally in Section |7| we determine the amount of energy radiated to infinity. 



3 Representation of the self-force 

In this section we show that the self-force F^(t) from ( |2.1| ) can be written 
in the form ( |2.2|) . We carry out this computation on the original fast time 
scale corresponding to ( |1 . 1| ) since we will need some of the arguments from 
0]. Thus we consider 

F s (t) = ( d 3 x(f)(x,t)Vp(x - q(t)). 
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Since <f>(x,t) = <p r (x,t) + (f> (x,t), where 0o = A0o with the initial values 
(po(x, 0) = and 7r (x, 0) = vr°(x), and since 

<pr(x,t) = [ 7^-/ d 2 yp(y-q(s)) 

AlT JO t — S J\y-x\=t-s 

is the retarded potential, we can decompose accordingly, 

F s {t) = F (t) + F r {t) = (0 O (-, t), Vp(- - q(t))) + t), Vp(- - q{t))} . 

Lemma 3.1 The function F (t) vanishes for t > t\ — 2R p /(l — v). 

Proof : Let U (t) denote the group generated by the free wave equation in 
L> 1,2 (IR 3 ) © L 2 (IR 3 ). Then (|D|) and Fourier transformation implies 

(<j>°(x), 7t°(x)) = - f° ds [U(-s)p(- -q°- v°s)](x) 



with p(x) = (0,p(x)). Thus Kirchhoff's formula yields, as a consequence of 
< 1, that 4>o(x, t) = for \x — q°\ < t — R p . Since \q{t) — g°| < vt, the 
claim follows. □ 

Hence to show ( fO| ) it is enough to prove 

Lemma 3.2 For t>ti, 

F r (t) = -m f {v(t))v{t) + a(v{t))v{t) + b{v{t),v{t)) + 0(e 3 ) , 

cf. fiT7]\), gJD, and Q). 

Proof: We follow the proof of || Lemma 5.1], but expand 

q(s) = q(t) - v(t)(t -s)+ l -v{t){t - sf - h(t)(t - s) 3 + 0{e 3 ) 



up to third order, which is allowed by Lemma |2Jj . Through Fourier trans- 
formation we arrive at 

FJt) = (-i) f ds ( d 3 k \p(k)\ 2 -^- sin |jfe|(t - s) e^^X*-*) 
Jo J \k\ 

xe -i\~\{k-v){t-s) 2 + \{k-v){t-sf] + ^ 
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with v = v(t), etc.. As in |5], Lemma 5.1], here and in the following / * ds(. . .) 
can be changed forth and back to //_ T ds(. . .) for all t, T > t\. Because 

e - ih i( fe . i ) (t - s )2 + i( fe .« )(t - s) 3 ] = 1 + i (A; . ^(t-sf-Uk-vMt-sf 

2 6 



for t — s = 0(1) by ( fg.l| ) below, we obtain, for t,T > ti, 

F r (t) = (-i) / T / ci 3 /Hp(fc)| 2 -^ sin|A;|Te-^ r 
Jo J |fc| 

1 + ~{k ■ v)r 2 - Uk ■ v)t 3 - hk ■ v) 2 r 4 } + 0(e 3 ) . 
2 o o J 

Let 

f T d sm\k\r i{k . v)T 
h J o ar e r , p U,...,4. 

Then 

(v ■ V v ) V v h = -k(k ■ v)I 3 and (v ■ V v ) 2 V v h = ik(k ■ v) 2 I 4 . 

Our claim now follows from Lemma [D| below, J <i 3 £; |/3(&;)| 2 &;/o — > 0, and 
(1/2) JcZ 3 A; |p(A;)| 2 A;(A; • v)I 2 -> -m f (v(t))v(t) for T -> oo; see f| Appendix 
A]. □ 

Lemma 3.3 VFe /iai>e t/ie identity 

' dtt [ d^kimf^^e-^ = (e 2 /4vr) 7 2 . 
J \k\ 

Proof: Since p(k) = p r {\k\) is radial, and by transformation to polar coor- 
dinates, 

,sinjfc]t i(k . v)t = 4?r /•«> 

1*1 ~ *h 

Thus for fixed T > 0, 



/ d 3 k\p{k)f-^^ e~ i{k - v)t = ^- f°° dR\p r {R)\ 2 sin{Rt) sm{Rt\v\) 

J \nj\ T \ V \ J 



3{.UfW|2^p-i(Hi 



t / effc|p(fc)| 
J \k\ 

2tt /•<*> di? 2 

= u / -p brC-R) 
f JO -ft 





u|)r) sin(i?(l + 


«|)T) 


11 


1 + 
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To complete the proof we only need to verify that / d 3 k\p(k) | 2 \k\ 3 sin \k\T — > 
e 2 /47r as T — > oo. To see this, let i[)(k) = \k\~ 3 sm(\k\T) . Then 

/ d 3 k\p{k)\ 2 ^j{k) = (2tt)~ 3/2 J d 3 xp{x) J d 3 y p{y)^{x - y) 

and we are going to show ip(x) — > \Jn/2 as T — > oo. We have, by transfor- 
mation to polar coordinates, 

(2tt) 3 /^(x) = [d 3 k^k)e- ikx = ^ r ds M^ Ms\x\m ^ 2 ^ 2 

J Jo s s\x\/T 

for T — > oo. This completes the proof. □ 



4 More about the center manifold 

In this section we explain the behaviour of solutions on the center manifold. 
First we show that the unperturbed manifold Iq from ( |2.6| ) is hyperbolic in 
normal direction. 

Lemma 4.1 The eigenvalues of D y g(x,y,0) = a{x2)~ 1 m(x2) are bounded 
below by a positive constant, uniformly in x = (r,u) with r E H 3 and \u\ < 
1 — 5, for all prescribed 6 e]0, 1]. 

Proof: By Thm. 2, p. 185], a(u) and m{u) can be simultaneously trans- 
formed to diagonal form through a single non-singular matrix B. In ad- 
dition, denoting by bj ^ the jth column of B and by Xj the jth eigen- 
value of a{u)~ l m{u), one has Xja(u)bj = m{u)bj, j = 1,2,3. Multiplication 
by bj leads to \ j (e 2 /12ir)j 3 [jb 2 j + 4j 3 (v ■ bj) 2 ] > 7& 2 + 7 3 (> ■ bj) 2 , and thus 
\j > (37r/e 2 ) 7 - 3 . ' " ' □ 

Since a{u), m{u) are modified to be constant outside \u\ < 1 — 5, their 
corresponding eigenvalues are uniformly bounded below for \u\ < 1. As a 



consequence of Lemma [O the manifolds X £ are unstable at some exponential 



rate e Mi for solutions in the normal direction. 
We note that, by JTTJ Thm. 2.1], 

sup{\h £ (r,u)\ : (r,u) etfxV.e e]0,e ]} <c = c(5). (4.1) 
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Our next aim is to prove global existence of solutions to ( |1.19| ) forward in 
time which start over /Q = IR 3 x {u G IR 3 : \u\ < v} on the center manifold, 
provided e < E\ with e\ > sufficiently small. For this purpose we introduce 
a suitable Lyapunov function. 

Lemma 4.2 Let 

G e (r, u, u) = H(r, u) — e(a{u)u) ■ u = £ s {u) + V(r) — e(a{u)u) ■ u . 
Then along solutions (r(t),u(t),u(t)) of ( \1.19j ) we have 

jG £ (r, u, u) = -e(e 2 /127r) [6 7 8 (u • uf + 7 6 u 2 ]. (4.2) 

Proof : Observing that 

(a(u)u) ■ u = (e 2 /127r) 7 6 (1 + 3u 2 )(u ■ u), 

this is a straightforward calculation. □ 

Through the Lyapunov function G e we can control the long time be- 
haviour. 

Theorem 4.3 Let (U) or (U') hold and let any global solution (r(t),u(t)) of 
( \1.13j ) be given such that sup t>0 \u(t)\ < u(e) < 1 and sup t>0 \it(t)\ < c(e), 
for possibly e-dependent constants u(e) and c(e). Then 

ii{t) 0, u{t) -> 0, and W(r(t)) ^0 as t -> oo. 

Proof : Denoting by c(e) or C (e) general e-dependent constants, by Lemma 
4~2] we have along a trajectory 

c(e) [ T u 2 dt < - [ T ^-G £ dt 
Jo Jo dt 

= -£ s (u(T)) - V{r{T)) + e(a(u(T))u(T)) ■ u(T) 

+£ s (u°) + V(r°) - e(a(u°)u°) ■ u° 
< C(e, data). 

For the last estimate observe inf rg]R 3 V(r) > — oo in both cases (U) and (V). 
Thus f °°u 2 dt < C(e, data) and, by ( |1.19| ), also sup t>0 \u(t)\ < C(e, data). 
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Hence we conclude it(t) — > as t — > oo. Next, differentiation of ( |1.19| ) yields 
sup t>0 | ii (t)\ < C(e, data), and thus from u(t) — > we find u(t) — > 0. 
Therefore Vy(r(t)) — > follows from the equation ( |1.19| ). □ 

In the demonstration of the following theorem we use the sublevel sets 
{G e < c} = {(r,u,u) : G e (r,u,u) < c} and {H < c} = {(r, u) : H(r,u) < 
c} for c G H. However, before proceeding, we first have to introduce an 
appropriate 5 = 5(v) > small to modify the functions a(u), m(u), and 
b(u, it) outside \u\ < 1 — 5, cf. Section |[ To do this, we assume (U) from now 
on. The case (If) is discussed in the remarks below. Since V is bounded and 
v < 1, we can find c G IR such that K,y C {H < c }. Then as a consequence 
of S s (u) — > oo for |tt| — > 1, we have 

So = sup : (r,u) G {-ff < Cq + 1} for some r G IR 3 | < 1 . (4.3) 
Let us define 

5 = min{(l-w)/2, (1 - s )/2} > 0. (4.4) 

Theorem 4.4 ^sswrne ^/ie potential V to satisfy the condition (U). Then 
there exists E\ > depending only upon v such that for e g]0,£i] all solutions 
of ( \l.iyj ) starting at points (r,u, h £ (r,u)) G X £ , \u\ < v, stay away from the 
boundary {(r,u, h £ (r,u)) : \u\ = 1 — 5} for all future times. In particular, 
solutions exist globally. 

Proof : Let us denote the bound c(S) from ( |4.1| ) by c\ and let us fix c a > 
such that \a(u)\ < c a for all \u\ < 1. We recall that a{u) was modified to be 
constant outside \u\ <1 — 5. We define Si = mm{e , (2c a Ci) -1 } > 0. 

Let (r, u) G /C s . Then G £ (r,u, h e (r,u)) = H(r,u) — e(a(u)h e (r,u)) -u < 
co+c a ciS. Because of Lemma the set {G e < co+c a cie} is forward invariant 
and the solution remains in this set for all future times. On the other hand, 
since v < 1 — 2<5 < 1 — 5, the solution of the modified problem is a solution 
to ( |1.19| ) and stays on X e , at least for a short times. For the fixed time span 



where this holds the solution is of the form (rj, U\, h e (r\, U\)) and we have 
F(ri,«i) = G e (r 1 ,u 1 ,h e (r 1 ,u 1 ))+e(a(u 1 )h e (r 1 ,u 1 ))-u 1 < c + c a ci£ + c a ci£ = 
Co + 2c a ciE < Co + 1 for e < E\. Therefore by ( [4.3|) , \u\\ < s < 1 — 25 < 1 — 5. 
This argument shows that in fact the solution is confined to {(r, u,u) : \u\ < 
1 — 25}. Hence the solution of the modified problem exists, is a solution to 
( |1.19| ), and stays on I e for all future times. □ 
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Corollary 4.5 In the setting of Theorem \!^.J\ , for solutions of starting 
on Z £ , 

sup{|u(t)| : t e H,e e]0,ei]} < 1 - 28 < 1 , and 
sup{|«(t)| :t G M,e e]0,ei]} + sup{|ii(t)| : t G M,e g]0,£i]} < c(5). 

(4.5) 



In particular by Theorem ^73 

ii(t) -> 0, -> 0, and W(r(t)) ^0 as t -> oo. 

Proof : The first estimate was mentioned already in the preceding proof. 
For the second we note that ( |4.1|) applies, since the trajectory stays on the 
center manifold, u(t) = h £ (r(t),u(t)). Concerning the last bound, we may 
write 

h £ (r, u) = — m(M) -1 W(r) + hx t£ (r, u) with \hi >£ (r,u)\ < c(5)e (4.6) 

for (r,u) G IR 3 x V, see JTT], Thm. 2.9]. By ( [TT9D , 

\eu\ < \a(u)~ 1 \\m(u)hi i£ (r,u) — eb(u,ii)\ < c(5)e, 

so we are done. □ 

Solutions on X £ are uniformly bounded, in the sense of the corollary; in 
general a bound on r(t) cannot be expected, e.g. in a scattering situation. 
Conversely, as to be shown next, solutions with uniformly bounded u(t), u(t), 
and u(t) are confined to the center manifolds. 

Proposition 4.6 Suppose we have a family (r £ (t) , u £ (t)) , e g]0, £2]? of solu- 
tions to ( \l.ldj ) such that 

sup{\u e (t)\ : t G IR,£ G]0,£ 2 ]} < u < 1 , and 
sup{|w £ (t)| : t G H,£ G]0,£ 2 ]} + sup{|ii £ (t)| : t G H,£ G]0,£ 2 ]} < c 2 . 

Then for sufficiently small e the solutions have to lie on T £ . 
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Proof: Note that we can construct X £ here by modifying a(u),m(u), and 
b(u,u) to be constant outside, say, {u : \u\ < (1 + u)/2}. According to 
Thm. 2.1 (ii)] there exists 5 > such that for all e small and solutions 
(x(t),y(t)) to (|2~5|) the condition sup t6R |w(t) — /i(x(t))| < 5 implies that the 
solution has to lie on X e . With x(t) = {r £ (t) , u £ (t)) and y(/f:) = u £ (t), this 
condition is verified since we obtain from ( |1.19| ) and the assumed bounds 
\u £ (t) + m{u £ {t))- l WV(r £ {t))\ <ce<5, the latter for e small. □ 

The asymptotic condition, r(t) — > 0, of Dirac and Haag is also sufficient 
for a solution to lie on I £ , in the following sense. 

Proposition 4.7 Suppose a family (r £ (t) , u £ (£)) , £ g]0, e 2 ] ; of solutions to 
( \1.19j ) is given such that 



sup{> e (t)| : t e E,£G]0,£ 2 ]} < u < 1 

and r £ (t) = u £ {t) — > as t — > oo for each e g]0, e 2 ]. T/ien /or sufficient small 
e the solutions have to lie on X £ . 

Proof: Fix 8 > 0. Since Theorem fO applies, we find in the notation of 



Proposition 4.6 



\y(t) - h(x(t))\ = \u £ (t) + m(u £ (t))- l VV(r £ (t))\ < 5/2 

for t > t(e), with some t(e). Thus the solution remains (5/2)-close to X 
after time t(e), and hence by ( |4.6|) also 5-close to X e for e small. Since X e 
is normally hyperbolic (repulsive) at an e-independent rate and since 5 > 
was arbitrary, this can only happen if the solution was already contained in 
X. ' " □ 



Corollary [4.5| provides a partial information on the long time behavior of 
the solutions to ( |2.5| ) on the center manifold. Roughly one can distinguish two 
classes, (i) (scattering): The particle enters a domain where — VV = at r*i 
with velocity Uoo. If the straight line trajectory r\ +Uoot, t > 0, is contained 
in this domain, then the particle travels freely to infinity. Physically this 
is a scattering trajectory. In this case lim^ 00 M(t) = ^ 0, whereas the 
position has no limit, (ii) (bounded motion): We assume that \r(t)\ < const. 
and that within this ball the critical points of V form a discrete set. Then by 
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Corollary [15] and by continuity we have lim^ 00 -u(t) = and lim^oo r(t) = 
Too, where is one of the critical points of V. If is a stable critical point, 
then the relaxation is exponentially fast, as can be seen from linearization 
around the fixed point. Clearly (i) and (ii) do not exhaust all possibilities. 
The critical points of V could lie on a sphere. If V is confining, one would 
still expect convergence to a definite r^. Moreover, V could vanish inside a 
ball. If — W is pointing towards the ball, then close to each turning point 
the particle looses energy. Thus lim^ 00 M(t) = 0, whereas the position has 
no limit. The potential could decrease so slowly at infinity that no definite 
velocity is approached. All these cases have to be studied separately. 

Up to now we discussed bounded potentials satisfying (U). In the intro- 
duction we claimed that our results remain valid also for confining poten- 
tials satisfying (U 1 ). In this case, since V is unbounded, we have no longer 
fCv C {H < Cq} for some Cq € E as in Theorem [4.4| above. However, by en- 



ergy conservation, one can derive the a priori bound sup 4g]R \q £ (t)\ < M for 
solutions to the true system ( |1.5|) on the macroscopic time scale. Thus the 
motion is bounded also in the g-direction and it suffices to build the center 
manifold for the effective equation (|1.19|) over the bounded domain 

/C^_ = {r G H 3 : |r| < M} x {u G H 3 : |u| < v} , 

enlarged to a suitable K<m+ii-5 such that solutions starting over JC^^ stay 
away from the boundary of Km+x 1-5 for £ > sufficiently small. In this 
manner we obtain 

Theorem 4.8 Assume (If) holds for the potential, and let K,^ g be defined 
as above. Then there exists E\ > depending only on the initial data such 
that for e g]0, ex] all solutions of starting at points (r,u, h £ (r,u)) G 

T £ , (r, it) G Km i, exist globally. Moreover, these solutions are uniformly 
bounded, 

sup{|r(t)| : t G H,e G]0,ex]} < c(S), 
sup{|«(t)| : t G IR,e e]0,e a ]} < 1 - 25 < 1, and 
sup{|u(*)| :t G M,e g]0, ei]} + sup{|ii(t)| : t G IR,e G]0,e x ]} < c{8). 

(4.7) 

In addition, 

u{t) -> 0, u{t) -> 0, and W(r(*)) ^0 as t -> oo. 
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Proof : The proof is similar to the one of Theorem [4.4| . Concerning the 
boundedness, note that again for some cq G IR and e > small, 

K-m,s C {H < c } C {(r, u) : G e (r, w, ft, e (r, u)) < c + c a cie} C < c + 1}. 

Thus all solutions starting over K^s wm remain on the manifolds over {H < 
Co + 1}. Since this set is independent of e and compact by (U'), the solutions 
must be uniformly bounded, because h £ is uniformly bounded. □ 

On the center manifold the motion is governed by the (second order) 
equation 

x = f(x,h £ (x)). (4.8) 

Since the existence of h £ is established only abstractly, Eq. (|4.8|) is somewhat 
implicit. From [jTT], (2.9-1) & Thm. 2.9] we know that h £ depends smoothly on 
e. Thus (|4.8| ) can be expanded in e. Including the first Taylor term we pick 
up an error of order e 2 , which is of the same order as the error between the 
true and the comparison dynamics on the center manifold. For consistency 
we should stop then at this order. We make the ansatz 

h £ (r,u) = h (r,u) + e/ii(r,u) + h 2 , £ {r,u) , \h 2 , e {r, u)\ < c(5)e 2 

for (r, u) G IR 3 x V. Then 

m(u)h (r,u) = — VV(r), 

and h\(r,u) is determined through 

D x h (r, u)f(r, u, h (r, u)) = D y g(r, u, h (r, u), 0)/ii(r, u)+D £ g(r, u, h(r, u),0), 

see |TT|, (2.9-1) & Thm. 2.9]. Computing the respective derivatives one arrives 
at 



m(u)hx(r, u) 



a[u) 



— m[u) 



_1 VV(r)u 



+ 



( 



-m[u 



Vl / (r),m(n)" 1 Vl / (r 



\du 

+ b{u,m(uy 1 WV{r)) 

and the effective second order equation 

r = u, m{u)u = —W(r) + sm(u)hx(r,u) 
of the particle motion on the center manifold. 



(4.9) 
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5 Comparison of the true and the effective 
system 



In this section we prove Theorem |2.2| . Since we have = | ^ £ (^1) | < v by 



( |1.8|) , Theorem |4.4| , resp. Theorem |4.8| , implies that the solution trajectory 
of the system with the modified functions a(u), m(u), and b(u,u) is indeed 
a solution trajectory to ( P-.19Q . Recall that, by (|1.18|) and ( |1.19|) , 

m(v £ )v £ = -W(g e ) + ea(v £ )v £ + eb(v £ ,v £ ) + e 2 f £ (t), t>et u (5.1) 
- W(r) + ea{u)u + eb(u, u), (5.2) 



m(u )u 



with |/ e (£)| < C. Using the bounds ( |TT8] ) and (fO>D, resp. (JO]), we infer the 
weaker estimate 



m(v £ )v £ 



m[u)u 



which has been proved already in ||. Hence 
there. 



\7V(q £ ) + <D(e), t>t u 
W(r) + 0(e), 

follows by the argument 



To show ( |2.8j ) we compute as in Lemma |47|, using (I57T 
d 



dt 



G £ (q £ (t),v £ (t),v £ (t)) 



e 2 f(t)v £ (t) 

-e(e 2 /12n) \j{v £ {t)fv £ {tf + 6~?(v £ (t)) s (v £ (t) ■ v £ {t)f 



t > eh. 



Since r(eti) = q £ (eti) and u(eti) = v £ (eti), using the uniform bounds, we 
have for t > et\ 

\H(q £ (t),v £ (t))-H(r(t),u(t))\ 



< e 



a{v £ {t))v £ {t)) ■ v £ {t) - (a(u(t))u(t)) ■ u{t) 



- I ds 

/etl 



e 2 \f £ (s)v £ (s)\+e(e 2 /12n)[ 7 ({; £ (s)) e t; £ (s) 2 - 7 (M(s)) e M £ (s) 



+6 



7 (v £ (s)r(v £ (s) ■ v £ (s)Y - j(u(s)r(u(s) ■ u(s)Y 



< Ce 



\v £ {t) -u{t)\ + \v £ {t) -u(t)\ +CeH 



v £ (s) - u(s)\ + \v £ (s) - ii(s)\ 



+Ce / ds 

Jeti 

< Ce 2 (l + t) <Ce 2 , 
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by ( p.7| ) for t = 0(1). This concludes the proof of Theorem |2T2 . 



□ 



Finally we show that on a microscopic time scale our results track the 
true trajectory with a higher precision, cf. (|2.9|), Remark |2~3"|(i). 

Proposition 5.1 We have 

\q £ {t) -r{t)\<ce 3 and \v £ (t) - u(t)\ < ce 2 , t = 0(s), 

i.e., holds. 

Proof: Define = (e~ 1 q £ (es) — e~ 1 r(es), v £ (es) — u(es), ev £ (es) — eu(es 

( 1 

for s > ti. Then *(s) = AV(s) + 6(s), where A = 1 | and 6(s) 



£ 2 (0, 0, v £ (es) — u(es)), whence 

|0(s)| < cs (\q £ (es) -r(ss)\ + \v £ (es) 



1 
V 



— u(ss)\ + \v (es) — u(es) 



< c(\^(s)\+e 3 ), s>t 



i- 



by (|5.1|) , (|5.2|) , and the uniform bounds. Therefore by the variation of 
constants formula and Gronwall's inequality for s G [ti, ti + r], |^(s)| < 
c(r)(|\E'(si)| +e 3 ). Consequently, q £ (sti) = r(et 1 ) and v £ (et 1 ) = u(eb\) yields 



e' l \q £ {t) -r(t)\ + \v £ (t) - u(t)\ < c(r) [e\v e (et x ) -«(e*i)| +£ 



for t e [eti,eti + er]. By Q and (|5"^) ; |w £ (rfi) - u(eti)| < ce, so that ( gjg ) 
follows. □ 



6 Adiabatic limit of the fields 



We prove Theorem 2A . Let U it) again denote the fundamental solution of 
the wave equation in D 1 ' 2 ^ 3 ) © L 2 (IR 3 ). We set Z{x,t) = {(j){x,t),-K{x,t)) 
as well as p £ = (0,p e ). Then the inhomogeneous wave equation in ( p..5| ) is 
solved as 

Z(x,t) = [U(t)Z(-,0)}(x) -y/i fds [U(t - s)p £ (- - q £ (s))}(x). 

JO 
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Since 



Z(x,0) = - v / i [° ds[U(s)p £ (--q -v°s)](x), 

J — oo 



cf. Lemma |3.1| , we have for t > 

Z(x,t) = -Vi f ds[U(t-s)p e (-~q e (s))}(x), 

J — oo 

where we extended the position to negative times t < by g £ (t) = q° + t>°£. 
Thus by the solution formula for the wave equation 

-L<f>( X ,t) = - f p- p £ (y- q e( t -\x-y\)) (6.1) 

and ii(x,t) = <fi(x,t). For £ — > 0, g £ (t) — > r(t), cf. (|1.9f) , with r(t) extended 
to negative times by r(t) = q° + v°t. Moreover, p e {x) = e^ 3 p(e^ 1 x) — > e^o 
in the sense of distributions. Hence the transformation z = y — q £ (t — \x — 
y\), det(dy/dz) = [1 — v £ (t — \x — y|) • (x — — y\]~ , in ( |6.1| ) yields the 
pointwise convergence ( j2.10|) , except on the worldline of the particle, since 
the integrand in ( |6.1| ) is singular at y — x, i.e. for x = r(t) which corresponds 

to £ re t = t. 

The analogous argument works for ir(x,t). In the limit e — > 0, tt is dis- 
continuous at the light cone {x : \x\ = t}, which we avoided due to our 
assumption. □ 



7 Radiated Energy 

Let E Rtq sM(t + R) be the energy, particle plus field, at time t + R in a ball 
of radius R centered at q e {t). For R > eR p this energy changes as 

-j t ^E R>q e {t) (t + R) 

= ^(n mac (t = 0)-^ { d z x\\n(x,t + R)\ 2 + \V<P(x,t + R)\ 2 \) 

at \ 2 J{\x~q£(t)\>R} 1 V 

= R 2 [ d 2 u7r(q e (t) + Ru,t + R) to ■ V0(<f(t) + Ru,t + R) 

J\u\=l 
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+- 



R 2 



w=l 



du (u ■ v £ (t)) [\7i(q £ (t) + Ruj, t + R)\ 2 

+ |V0(g e (t) +Rw,t + R)\' 



(7.1) 



where we used that the total energy is conserved. 

Er changes because there is energy flowing back and forth between par- 
ticle and field, and because energy is lost irreversibly to infinity. To separate 
both contributions we take the limit R — > oo. Using ( |6.1|) and the relation 
t + R- \q £ (t) + Ruj - y\ = t + uj ■ (y - q 6 (t)) + 0(1/ R) for bounded \y\, we 



arrive at 



I s it) 



lim ^-\E 



R-^oo dt 



(t + R) 



M=i 



d 2 u(l -uo-v £ (t)) 



x 



d 3 yp e (y-q £ (t + u;- [y - q £ (t)})) 
u-v £ (t + uj- [y -q £ (t)}) 



cf. |7|, Sec. 3] for details on a similar calculation. In fact, there the ball of 
radius R was centered at the origin and the second summand in ( |7. 1| ) is 
absent. To let e — > 0, we again transform to z — y — q £ (t + uj ■ [y — q £ (t)}), 
det(dy/dz) = [1 — uj ■ v £ (t + uj ■ [y — g e (t)])] _1 , use p £ (x) — > e^o in the sense 
of distributions, and insert the identity y = q £ (t) for z = to obtain 



lime" 1 I £ (t) 



-e 2 (4n)- 2 [ d 2 u (1 - uo ■ u(t))- 5 (uj ■ ii(t)f 

J\u\=l 

-(e 2 /12n) [6 7 8 (n(t)-«(t)) 2 + 7 6 «(t) 2 ], 



in agreement with 

Alternatively, we could first take the limit e — > in ( |7.1| ). Using Theorem 
2.4 we find, with (0, w) denoting the limit fields from ( |2.10|) , (|2.11|) , 



Hi) 



lim s 1 ^- ( 

e^O dt V 



E 



(t + R) 



Rf 



|«|=1 



d 2 uj W(r(t) + Ruj,t + R)uj ■ V<p(r(t) + Ru,t + R) 



R 2 

+ — / d 2 uj (uj ■ u(t)) \T(r(t) + Ru } t + R)\' 



UJ =1 



+\V(f>(r(t) + Ruj,t + R)\' 



25 



Since both W and V0 have one term proportional to Br 1 and other contribu- 
tions of order R~ 2 , in the limit R — > oo only the product of the two leading 
terms survives and it follows that 

lim I R (t) = -e 2 (4vr)- 2 / d 2 u (1 - oo ■ u(t))- 5 (lu ■ ii(t)) 2 , 

R^oo J\u\=l 

as before. We note that the radiated energy is of order e and it therefore 
suffices to use the effective dynamics to order one, i.e., ignoring the radiation 
reaction. 



8 Appendix: Proof of Lemma gTT] 

In this appendix we prove Lemma |2.1|. Since we need to use some identities 



from ||, we switch back to the original time scale of Hence we have 

to show 



Lemma 2.1 For solutions of (1.1 ) with initial values satisfying (1-4), i-e., 
starting on the soliton manifold, and for \p\ L 2 sufficiently small we have 

sup | v (t)\ < Ce 3 . 
ten 

The constant C and the bound on \p\ L 2 depend only on the data. 

Proof: From [^|, Lemma 2.2 and Prop. 4.1] we already know the bounds 

suPteR \v(t)\ <v<l, sup teR \v(t) \ + sup 4eR \p(t)\ < Ce 

and sup teR \v(t) \ + sup te]R \p(t)\ < Ce 2 , (8.1) 



for \p\ L 2 sufficiently small. The constants v and C appearing in ( |8.1| ) do not 
depend on the particular solution, but only on bounds for the initial values. 
Denote 



<p(x,t) \ _ / <j>{x,t) - <f> v ( t )(x - q{t)) 
ip(x,t) J \ir(x,t)-ir v (t)(x-q(t)) 



Z(x,t) = 

Then, cf. §, with L{t)<f> = V0 • v(t) + <p, 

p(t) = -e 2 V 2 V{eq{t)) ■ v(t) + J d 3 x (L{t)<p)(x + q(t),t)Vp(x) 
=: -e 2 V 2 V(eq(t)) ■ v(t) + M(t). 
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Therefore 

p (t) = - e 2 V 2 V(eq(t)) ■ v(t) - e 3 V 3 V(eq(t))(v(t),v(t)) + M(t) . (8.2) 
Below we will show 
Lemma 8.1 The estimate 

|MWI<c( £ 3 + |p| L! { T I|ML A ) 

holds. 

Then according to (|8.2|) , (|8.1| ), and assumption {U) on the potential, 

.3 . . I* \v (s) 



Since 



< C(\'P\+e 3 ), 

the claim of Lemma 2A obtains from ( |S.3|) by taking \p\ L 2 small enough. □ 



nt)\<C[e> + \p\„] a Y^^ds). (8.3) 



Thus it remains to give the 
Proof of Lemma |8.1j : First note 

M(t)= ( d 3 x((£(t)Z(;t))(x),Vp*(x-q(t))) p*{x) = {p{x),0), 



where C(t)Z = VZ ■ v(t) + (V 2 Z)(v(t), v(t)) + 2VZ ■ v{t) + Z. Because 
Z = AZ - B, with 

a„d «m).(^<::j|;»:|>), M 

we obtain 

j t (£(t)Z) = A(C(t)Z) - C(t)B + 2 [(V 2 Z)(w, u) + VZ ■ v] + VZ • i) . 
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Let U(t) again denote the group generated by the free wave equation on 
_D 1,2 (IR 3 ) © L 2 (IR 3 ). Then 

M(t) = (U(t)[C(0)Z(.,0)],Vp^-q(t))) LHR2) 



+ / ds 
Jo 



-(U(t-s)[jC(s)B(.,s)],Vp*(--q(t)) /L2ii2) 



+ 2 (U(t - s)[(\7 2 Z(; s))(v(s),v(s)) + VZ{; s) ■ v(s)}, 



'L 2 (IR 2 ) 

+ (U(t - s)[VZ{; s) ■ v{s)\, Vp*(- - q(t)) /L2{n2) 
=: To + Tx + Ta + Ts. 

We estimate each term Tj separately, keeping all parts which do contain only 
initial values. Note that here according to (L4) we have Z(x, 0) = 0, so all 
these terms vanish. Nevertheless, we wanted to derive the general form of 
the estimate; see Remark |2.3| (iii). 

Estimate of T 3 : Since Z = AZ — B, we find 

Z(t) = U(t)Z(0) - /* ds U(t - s)B(-, s), (8.5) 
J 

and hence 

^ dsllKtMV Z (■.(]) -v(s)]. VnJ- -a(t))\ 

L 2 (IR 2 ) 



T 3 = J o ds(U(t)[VZ(.,0)-v(s)},Vp*(--q(t)) 

- J* ds J° dr (U(t - r)[VB(; r) ■ v(s)), Vp*(- - q(t)) 



T3.0 + T 3 i . 



L 2 (R 2 ) 



Then Lemma |8.2j below and (|8.1|) imply through integration by parts in the 
G? 3 x-integral, 

T 3A < Ce 2 f ds f dr < Ce 3 . 

Jo Jo 1 + (t — r) d 



Estimate of Tq: This term is determined solely through the data. 
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Estimate of T x : If we calculate the form of C(t)B = VB-v + (V 2 B)(v,v) + 
2VB ■ v + B explicitly from (|8.4j), fortunately many terms cancel, and we 



find with *.= ( 7 I, 

C{t)B = V v $ v (x - q)- v +3V 2 $„(x - q)(v,v)+Vl$ v (x - q)(v,v,v) ■ 



Now we may argue analogously to the estimate of T 3 and Lemma to 
obtain with 

( fa) ' 

V %(x) / 

the estimate 



[U(t-s)C(s)B(;s)}(x) 



\V(j>(x + q(t))\ < j 



C 



+ (t- s) 2 

Here we have used (|8.1|) and some of the estimates 



(e 3 + | v (s)|) , \x\<R p , t>s. (8.6) 



|VV^ 


+ |W 2 0, 


,{*) 


| + |VV 3 0„(x)| 


< 


C(l + 


a; )" 


-2 


\V 2 V v <f) v 


;x)| + |v 2 v^ 


»\ 


+ |V 2 V 3 0,(x)| 


< 


C(l + 


|x ) 


-3 


|V 3 V„0„ 


»| + |v 3 v^„ 




+ |V 3 V 3 t) (x)| 


< 


C(l + 


|x )" 


-4 


|V 4 V„i 


;x)| + |v 4 v 2 0„ 


»\ 


+ |V 4 V 3 ^(x)| 


< 


C(l + 




-5 


IVV^tt^j;)! + |VV 2 7T, 




| + \VV 3 v tt v (x)\ 


< 


C(l + 


\x\)' 


-3 


|V 2 V„7T, 


(*)| + |v 2 vX 


[x) 


+ |V 2 V 3 7T,(X)| 


< 


C(l + 


\x\y 


-4 


|V 3 V t) 7T t) 


(x)| + |V 3 VX 


[x) 


+ \V 3 V 3 v ir v (x)\ 


< 


C(l + 


\x\y 


-5 



for x G H 3 and |i>| < t>. From 



we conclude 



ds / d A x V(f)(x+q(t))p(x) < C e 3 +|p| 

J\x\<R 



L- 



V [S] 



o i + (t- s y 



1.7) 



ds | . 



Estimate of T 2 : Let P(t)Z = V 2 Z(-,t)v(t) + VZ(-,t). Then 
d 



dt 



(P(t)Z) = P(t)Z + (V 2 Z)v = A(P(t)Z) - P(t)B + (V 2 Z)t) . 
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Therefore by definition of T 2 , 

T 2 = 2 J*ds{U(t)[(P(0)Z(;0))-v(s)}, Vp,(--g(t))) i2(R2) 

+2 j'ds j* dr (U(t - r)[- P{t)B{-,t) + {V 2 Z(-,t))v{t) ■ r(.sj 
Vp*(--g(t)) 



'L 2 (R 2 ) 
= : ^2,0 + ^2,1 + ^22 ■ 

To estimate T 2>1 , observe 

p(t)S = VV^„(x - g) • v + VV^ 8 (i - g)(t>, v) . 

Hence we may argue as before to find |T2,i| < Ce 3 . In order to bound T 2i2 , 
similarly to the estimate of T 3 we again use (|8.5| ) to get 

^2,2 

= 2 jf*<fc J 3 dr(u(t)[V 2 Z(0)(v(r),v(s))},Wp^-q(t))) L2{u ^ 
r t 

-2 







ds J dr J d<7([/(t-<7)[V 2 5(-,<7)(v(T),i<s))], 



- 2~2 2.0 + T 2 . 2 . 



By ( |3.7| ) and the argument of Lemma [T2] then 

rt rs rr C £ 3 

T 2 , 2 ,i < ds dr da — — rr < Ce 3 . 

Jo Jo Jo 1 + (t — ap 

Summarizing all above estimates for T0-T3, we hence arrive at 

m)\ < c^+^jl^^^ds) 

+ {U(t)[C(0)Z(;0)}, Vp*(--q(t))) L2[n2) 

+ 2 jf* ds (U(t)[(P(0)Z(; 0)) ■ «(*)], Vp»(- " 9(*))> i2(]R2) 
t 



+ 2 /o ^ /o dl -(^(*)[V 2 ^(0)(v(r),ti(a))], Vp*(- - g(*))) L2(R2) 
+ jf (C/(t)[VZ(-, 0) • «(*)], Vp*(- - q(t))) . (8.8) 
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Concerning the terms that contain data, these vanish here since Z(x, 0) = 
as a consequence of (|1.4j) . This completes the proof of Lemma jSll . □ 



In case of solutions starting not on, but close, to the soliton manifold 
as discussed in Remark |273|(iii), conditions on the data have to be imposed 
to ensure the last four terms in (|8.8|) can also be estimated by Ce 3 . In |5|, 
Thm. 2.6] and Section 4 of that paper details are carried out for derivatives 
of one order less. 

Above we used the following lemma. 

Lemma 8.2 The estimate 

\\V[U(t - r)VB(-,r)](- + g(*))||^ < C 1 + , t>r, (8.9) 

holds. 

Proof: Such estimates have already been used in ||, but we nevertheless 
include some details of the argument. Let 



4>(x) 

7t(x) 



[U(t-r)VB(.,r)](x) 



for fixed t, r. By Kirchhoff 's formula for the solution to the wave equation 
and by O, 



d 2 y {t-T)V 2 V v n v{T) (y-q(r))-v(T) 



V(j)(x + q(t)) 
1 



4vr(t - r) 2 J\y~x-q(t)\ = (t-T) 

+V 2 V v <p v{T) {y-q(T))-v(T) 

+V 3 V ' v (f> v{T )(y - q(r))(v(r),y- x - q(t)) 

(8.10) 

Now |x| < R p and \y — x — q(t)\ = (t — r) yields \y — q(r)\ > (t — r) — v(t — 
T )-R p = (1 -v){t-T) -R p by (|8~T| ). As a consequence of (ET7D , hence (gg) 
follows from ( gig ). □ 

Acknowledgement: We thank A. Komech for useful discussions. 
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